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Abstract 



We report a measurement of the average _B-hadron lifetime using data col- 
lected with the SLD detector at the SLC in 1993. An inclusive analysis 
selected three-dimensional vertices with _B-hadron lifetime information in a 
sample of 50k Z° decays. A lifetime of 1.564 ± 0.030 (stat) ± 0.037 (syst) ps 
was extracted from the decay length distribution of these vertices using a 
binned maximum-likelihood method. 

Measurements of the i?-hadron lifetime tb are useful in exploring the physics of 6-quarks, 
particularly in determining the weak couplings of the b to lighter quarks. Precise measure- 
ments of the average value of tb remain interesting in view of the significant variation in the 
world average over the past few years (see Refs. [1-3]). Previous determinations of tb have 
relied on either tight lepton momentum cuts [4] or stringent vertex constraints [5] to isolate 
B hadron decay vertices. In the method presented here, all decay modes are used with high 
efficiency. 

This letter presents a measurement of tb based on inclusive three-dimensional recon- 
struction of secondary vertices in Z° — > bb events. The data for this analysis were collected 
at the SLAC Linear Collider (SLC) with the SLD experiment in 1993. The analysis uses a 
topological technique to select vertices with lifetime information and relies on Monte Carlo 
(MC) modeling to extract tb, using a maximum-likelihood method. This yielded a precise 
measurement of tb, even with a relatively small data sample. 

During the 1993 run, SLD recorded 1.8 pb~ x of e + e~ annihilation data at a center-of- 
mass energy of 91.26 ± 0.02 GeV. Charged particle tracking was provided by the central 
drift chamber (CDC) [6] and by the vertex detector (VXD) [7], with a combined impact 
parameter resolution in rcf>(rz) parameterized as a = 11(38) © 70/pV sin 6 /mi, where p is 
expressed in GeV/c. The liquid argon calorimeter [8] was used for triggering, thrust axis 
determination, and jet finding. The ( rms ) xyz profile of the SLC beams was approximately 
2.4 X 0.8 X 700 fim 3 at the interaction point (IP). The x and y positions of the IP were 
continuously measured, using reconstructed tracks from ~ 30 sequential hadronic Z° decays, 
giving a T X y = 7 ± 2 /mi [9]. The z position was measured on an event-by-event basis using 
the median z position of tracks at their point-of-closest-approach to the IP in the xy plane, 
with a resolution of 52 /mi for Z° — > bb events [9]. 

A detailed simulation of the detector and physics processes was used in this analysis. 
Hadronic Z° decays were generated using JETSET 6.3 [10] adjusted to reproduce data 
from other e + e~ experiments. The fragmentation function for b- and c-quarks followed 
the Peterson parameterization [11] with e& = 0.006 and e c = 0.06, respectively. A detailed 
description of the i?-hadron decay model may be found in Ref . [9] . Beam related backgrounds 
and detector noise were simulated by overlaying random trigger events that occurred in close 
time proximity to a Z° decay. Detector response was simulated with GEANT 3.15 [12]. 

Hadronic Z° decays were selected by requiring at least seven reconstructed tracks, total 
track energy greater than 18 GeV, and |cos Othrust | < 0.71 (thrust axis within CDC-VXD ac- 
ceptance). A sample of 29,400 events was selected with a nonhadronic background estimated 
to be < 0.1%. 
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A set of quality tracks for use in heavy quark tagging and vertexing was selected. Tracks 
measured in the CDC were required to have > 40 hits have a hit at a radius r < 39 cm, have 
a transverse momentum p xy > 0.4 GeV/c, extrapolate to the IP within 1 cm in xy and 1.5 cm 
in z, and have a good fit quality (% 2 jd.o.j '. < 5). Tracks were required to have at least one 
associated hit in the VXD, and a combined CDC/VXD fit with y^ 2 jd.o.j '. < 5. Tracks with 
a 2-D (r<f>) impact parameter 6 > 3 mm or a 2-D impact parameter error as > 250 fim were 
removed. Tracks from identified 7 conversions, K° or A decays, were also removed. A dis- 
crepancy in the fraction of tracks passing quality track selection between data and MC was 
corrected by applying a momentum and angle dependent correction of ~ 6% to the MC [9]. 

A sample of 4299 events with an enriched Z° — > bb content was tagged by choosing 
hadronic events that have three or more quality tracks with normalized 2-D impact pa- 
rameter S/(Ts > 3.0. From the MC, this algorithm was determined to be 60% efficient at 
identifying Z° — > bb events providing a 90% pure sample for a sample with tb = 1.510 ps. 
This tag had a minimal effect on the measured lifetime due to a subsequent cut on flight 
distance, as will be discussed later. 

Secondary vertex reconstruction and selection proceeded in three stages. In the first 
stage, quality tracks with p > 1.0 GeV/c were used to construct secondary vertices. Two- 
prong secondary vertices were constructed from all pairs of these tracks in an event hemi- 
sphere, defined with respect to the thrust axis, if they extrapolated to within three standard 
deviations of a common point in that hemisphere. If a track was shared between a pair 
of two-prong vertices, a three-prong vertex was also constructed. Four-prong vertices were 
constructed in a similar way. Tracks from three- and four-prong vertices were required to 
extrapolate to within three standard deviations of a common point. The vertex fit proba- 
bility was required to be > 5%. To reduce the number of vertices containing tracks from 
the IP, the distance between each vertex and the IP was required to be greater than 1 mm, 
and at least one track in each vertex was required to have S/ag > 2.5. As a result, 97.8% of 
tagged events contained at least one selected vertex. 

On average, 6.5 vertices per hemisphere remained at this stage of the analysis. Figure 1(a) 
shows the number of vertices-per-event passing these selection criteria. The first column of 
Table I shows the composition of the vertices as determined from the MC. A vertex was 
classified as 'b ' if all the tracks came from the weak-decay vertex of a B hadron, and as a 
'b + cascade c ' if at least one track came from the weak decay of a B hadron and the rest 
from the weak decay of a hadron containing a cascade c-quark. The other categories were 
defined similarly. 

Because of the loose vertexing criteria, most geometric vertices were reconstructed; how- 
ever, tracks were allowed to contribute to more than one vertex. In the second stage of 
vertex selection, all possible sets (partitions) of independent vertices in each event hemi- 
sphere were found; i.e., no vertices shared tracks. Events were rejected if the total number 
of partitions exceeded 1000. This removed 1.0% of tagged events. For each hemisphere, 
we selected the partition with maximum value of the product M = Y\aii vertices P(x 2 ^ d.o.f.), 
where P(x 2 ? d.o.f.) represents the vertex fit probability and the product includes all vertices 
in a given partition. The partition selection criterion was chosen to provide a high efficiency 
for finding vertices with lifetime information and low backgrounds. The true i?-hadron decay 
vertex was not necessarily reconstructed. Column 2 of Table I shows the vertex composition 
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Figure 1. Distribution of (a) number of vertices per event, (b) M, (c) number of prongs per 
vertex, (d) total vertex momentum for data (points) and MC (histogram) with t b = 1.510 ps, 
at different stages of vertex selection. 



TABLE I. Vertex type in sample according to Monte Carlo. 





Stage 1: after initial 


Stage 2: after 


Stage 3: 


Vertex type 


vertex selection (%) 


partition selection (%) 


final sample (%) 


b 


19 


21 


23 


Cascade c 


14 


24 


26 


b + cascade c 


48 


35 


37 


b + other 


9 


5 


3 


Cascade c + other 


5 


6 


3 


Primary c 


3 


5 


6 


IP 


0.5 


1.3 


0.4 


Other 


2 


3 


2 


Vertex multiplicity 
per hemisphere 


6.5 


0.7 


0.6 
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after selection of the best partition, and Figure 1(b) displays the distribution of the product 
M. To further reduce the background, particularly from the IP, the angle between the vertex 
line-of-flight and the nearest jet axis was required to be less than 150 mrad, and the vertex 
tracks were required to have transverse momentum with respect to the vertex line-of-flight 
greater than 0.07 GeV/c. A sample of 5856 vertices (0.69 vertices/hemisphere on average) 
remained at this stage of the analysis. Figure 1(c) shows the distribution of the number of 
prongs per vertex. 




Figure 2. Decay length distribution for vertices passing all selection criteria in data (points) 
and MC (histogram). The MC distribution corresponds to that with the best-fit lifetime. 

In the third stage of vertex selection, one vertex was selected in each hemisphere in 
order to avoid multiple counting. For hemispheres containing more than one vertex, the 
vertex closest to the IP was selected to enhance contributions from b vertices over cascade 
c and other vertices. As a result, 60% of hemispheres in the 4299 tagged events contained 
a selected vertex. Of these, 92% had a selected vertex containing tracks associated with 
the weak decay of the B hadron. Figure 1(d) displays the total vertex momentum 
distribution at this stage. Column 3 of Table I, shows the vertex composition in the final 
sample. Since this analysis is sensitive to the accuracy of the MC modeling, we have verified 
that data and MC matched closely at the various stages of vertex selection (see Fig. 1). 
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The i?-hadron lifetime was extracted from the decay length distribution, i.e., the distance 
between secondary vertices and the IP (see Fig. 2). A tb dependence was introduced into 
the MC by reweighting the decay length distribution of the i?-hadron component. The MC 
contained two sets of Z° — > bb decays. In the first set, B mesons (baryons), representing 
91.1% (8.9%) of all B hadrons, were generated with an average lifetime of 1.55 ps (1.10 ps), 
and the average lifetime for B mesons (baryons) was 2.00 ps (1.42 ps) in the second set. 
The average value of tb in the two MC sets was thus 1.510 ps and 1.948 ps, respectively. 
A binned likelihood was computed for average lifetime values ranging between 0.7 and 2.3 ps; 
bins in the decay length distribution were combined to have a minimum of ten entries. Using 
all MC samples, the maximum-likelihood method yielded an average i?-hadron lifetime of 
t b = 1.564 ± 0.030 ps with a x 2 /d.o.f. = 27.9/20. 

Systematic checks of the method and the detector modeling were performed. Independent 
samples of the tb = 1.510 ps MC were used as data yielding lifetimes in agreement with 
the generated lifetime. The lifetime was also extracted from the data using the tb = 1.510 
ps and tb = 1.948 ps MC Z° — > bb samples separately, yielding values consistent with the 
reported result. In the two MC sets, the average values of the decay length distributions for 
each vertex type at the same reweighted tb were found to agree to within 1%. Variations in 
the measured lifetime between different epochs of the run, i.e., under different run conditions, 
were found to be within the statistical error of our quoted value. No systematic dependence 
on azimuthal angle could be inferred from the lifetime measurements from four different 
azimuthal sections of the detector. We verified that the initial heavy-quark tag had negligible 
effect by comparing the vertex composition and decay length distributions for tagged and 
untagged Z° — > bb events that passed all other cuts in the analysis. The 1 mm decay length 
cut reduced the backgrounds, especially from Z° — > cc events, and was found to reduce the 
sensitivity to the choice of tagging criteria. The tag was also investigated by considering 
only hemispheres opposite jets which were tagged by requiring at least 3 tracks in the jet 
with S/(Ts > 3.0. This yielded an average lifetime of tb = 1.616 ± 0.052 ps, which was 
independent of the heavy-quark tag and consistent with the reported result. We verified 
that the choice of vertex in the last stage of the analysis had little effect on the lifetime 
measurement. 

The systematic error due to detector modeling (see Table II) was dominated by the 
modeling of the reconstruction of quality tracks in the detector. This effect included an 
overall normalization between data and MC after quality track selection, as well as a 0.3 track 
variation in quality track multiplicity between different run periods [9]. Track parameters 
were varied to reflect uncertainties in the radial and longitudinal alignment within the VXD 
and to account for the differences between data and MC in the tails of the IP position 
distribution along the z-axis. The IP x and y positions were Gaussian-smeared by 100 fim 
for 0.25% of the MC events to reproduce tails observed in the data. 

The systematic errors due to the heavy-quark physics modeling are also shown in Table II. 
Systematic errors due to b and c fragmentation were determined by varying e& and e c in the 
Peterson function used in our MC [11], according to (xe) = 0.700 ±0.011 for b fragmentation 
and (xe) = 0.49 ± 0.03 for c fragmentation [13]. The b- fragmentation systematic error also 
includes a contribution due to the uncertainty in the shape of the fragmentation function, 
which was estimated by using the Bowler-Lund parameterization [14]. Sensitivity to the 
charm content of i?-hadron decays was estimated by taking into account the uncertainty in 
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the B° and B + branching ratios [3]. The uncertainties were inflated by 100% to account 
for the fact that no data are available for B s and B baryon branching ratios. The absolute 
fraction of B baryons in Z° — > bb events was varied by ±0.05. The absolute fractions of D + 
and charm baryons in Z° — > cc events were varied by ±0.037 [15] and ±0.03 respectively. 
Rb and R c [16] were varied by the uncertainty in the world average [3]. The multiplicity of 
tracks from i?-hadron decays was varied by ±0.3 tracks to reflect the current uncertainty in 
this value [17]. The MC momentum spectra of D + and D° mesons from i?-hadron decays 
were adjusted to agree with recent CLEO data [18]. The lifetimes of D + , D° , D S} and A c 
were varied by the uncertainty in their world averages [3]. 



TABLE II. Summary of systematic errors. 



Error source 


At b (ps) 


Quality track multiplicity 


0.013 


VXD alignment 


0.005 


IP X y position tails 


<0.001 


b fragmentation 


0.018 


c fragmentation 


0.005 


Charm content of B decays 


0.014 


_B-baryon fraction 


0.004 


Charm-hadron fraction 


0.012 


R b 


<0.001 


R c 


0.002 


_B-hadron decay multiplicity 


0.009 


_D-momentum spectrum in B decays 


0.010 


Charm-hadron lifetime 


0.004 


Minimum decay length in lifetime measurement 


0.006 


Decay-length binning 


0.007 


Decay-length cut in vertex selection 


0.010 


Vertex and partition limits 


0.003 


Monte Carlo statistics 


0.008 


Total systematic errors 


0.037 



Decay-length bin size was varied from 0.25 mm to 2 mm, and the minimum decay length 
used to extract the lifetime was varied from 1 mm to 3 mm to determine the variation in the 
measured lifetime. Samples passing different decay-length cuts at the vertex reconstruction 
stage were fitted over the same range to assess the effect of this cut. The systematic error 
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due to the cut on the maximum number of partitions was estimated to be negligible. This 
was accomplished by weighing the i?-hadron decay-length distribution in the Monte Carlo 
final sample to reproduce the distribution prior to the cut. 

In conclusion, an inclusive vertexing technique was used to determine the average B- 
hadron lifetime from a sample of 50 k Z° decays. The measured lifetime was 1.564 ± 
0.030 (stat) ± 0.037 (syst) ps, consistent with recent results from other experiments [4,5]. 

We would like to thank the personnel of the SLAC accelerator department and the 
technical staffs of our collaborating institutions for their outstanding efforts. 
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